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Spectrum of boson fields and two-point correlators are analyzed in a quantum bar system (a su- 
perlattice formed by two crossed interacting arrays of quantum wires), with short range interwire 
interaction. The standard bosonization procedure is shown to be valid, within the two wave ap- 
proximation. The system behaves as a sliding Luttinger liquid in the vicinity of the V point, but its 
spectral and correlation characteristics have either f D or 2D nature depending on the direction of 
the wave vector in the rest of the Brillouin zone. Due to interwire interaction, unperturbed states, 
propagating along the two arrays of wires, are always mixed, and the transverse components of the 
correlation functions do not vanish. This mixing is especially strong around the diagonals of the 
Brillouin zone, where transverse correlators have the same order of magnitude as the longitudinal 
ones. 



I. INTRODUCTION AND OVERVIEW 

Diverse D — 1 dimensional objects embedded in D di- 
mensional structures were recently investigated experi- 
mentally and analyzed theoretically. Rubbers and var- 
ious percolation networks are examples of such disor- 
dered D — 1 objects, whereas self-organized stripes in 
oxicuprates, manganites, nanotube ropes and Quantum 
Hall systems are examples of periodic structures of this 
kind. In some cases, the effective dimensions of such or- 
dered structures may be intermediate e.g. between D = 1 
and D = 2. They are especially promising candidates for 
studying novel electronic correlation properties, which, in 
particular, are relevant for the search of Luttinger liquid 
(LL) finger-prints in two dimensions. This challenging 
idea is motivated by noticing some unusual properties of 
electrons in Cu-0 planes in High-T c materials^. How- 
ever, the Fermi liquid state seems to be rather robust in 
two dimensions. In this respect, a 2D system of weakly 
coupled ID quantum wiresau looks promising. Indeed, a 
theoretical analysis of stable LL phases was recently pre- 
sented forj-apSystem consisting of coupled parallel quan- 
tum wirestTQ and for 3D stacks of sheets of such wires 
in parallel and crossed orientations! . In most of these 
cases, the interaction between the parallel quantum wires 
is assumed to be perfect along the wireta, whereas the in- 
teraction between the modes generated in different wires 
depends only on the inter-wire distance. Along these 
lines, generalization of the LL theory for quasi 2D (and 
even 3-D) systems is reported in Ref. |s| where the interac- 
tion between two crossed arrays of parallel quantum wires 
forming some kind of a network, depends on the distance 
between points belonging to different arrays. As a result, 
the grid of crossed arrays retains its LL properties for the 
propagation along both subsets of parallel wires, whereas 
cross-correlations remained non-singular. This LL struc- 
ture can be interpreted as a quantunijanalog of a classical 
sliding phases of coupled XY chainsEI. A special case of 
2D grid where the crossed wires are coupled by tunneling 
interaction is considered in Refs. H|U1 



In the present paper, a different course is elaborated. 
We ask the question whether it is possible to encode 
both ID and 2D electron liquid regimes in the same sys- 
tem within the same energy scale. In order to unravel 
the pertinent physics we consider a grid with short-range 
inter-wire interaction. This approximation might look 
shaky if applied for crossed stripe arrays in the cuprates. 
On the othet, hand, it seems natural for 2D grids of 
nanotubesEllEj, or artificially fabricated bars of quantum 
wires with grid periods a± t 2 which exceed the lattice spac- 
ing of a single wire or the diameter of a nanotube. It will 
be demonstrated below that the short-range interaction 
with radius ro <C ai t 2 turns out to be effectively weak. 
Therefore, such a quantum bar (QB) retains the II? LL 
character for the motion along the wires similarly to the 
case considered in Ref. ||. At the same time, however, 
the boson mode propagation along the diagonals of the 
grid is also feasible. This process is essentially a two- 
dimensional one, as well as the shape of the Brillouin 
zone and equipotential surfaces in the reciprocal QB lat- 
tice. 

Before developing the formalism, a few words about 
the main assumptions are in order. Our attention here 
is mainly focused on charge modes, so it is assumed 
that there is a gap for spin excitations. Next, we are 
mainly interested in electronic properties of QB which 
are not related to simple charge instabilities like com- 
mensurate CDW, so that the (for simplicity equal) pe- 
riods a\ — 02 — a are supposed to be incommensurate 
with the lattice spacing. The Brillouin zone (BZ) of the 
QB supcrlattice is two-dimensional, and the nature of ex- 
citations propagating in this BZ is determined by Bragg 
interference of modes with the superlattice wave vector. 
This interference (Umklapp processes) is, of course de- 
structive for LL excitations with both wave vector com- 
ponents close to multiple integers of 2ir/a. However, in 
case of weak scattering Vag/TiVF 1, where V is the in- 
terwire interaction strength, g is the usual LL parameter 
and vf is the Fermi velocity, only two- wave interference 
processes near the boundaries of the BZ are significant. 



f 



One can then hope that the harmonic boson modes sur- 
vive in the major part of the BZ, and that the Hamilto- 
nian of the QB might still be diagonalized without losing 
the main characteristic features of the LL physics. 



II. FORMALISM: HAMILTONIAN AND MAIN 
APPROXIMATIONS 

We consider a 2D periodic grid consisting of two arrays 
of ID quantum wires of length L oriented along unit 
vectors ei.2 with an angle if between them. The full 
Hamiltonian of the system is, 



H — Hi + H 2 + Hi, 



(1) 



The Hamiltonian Hi describes the ID boson field in the 
«-th array of wires (i = 1,2) 
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Here Vq is the coupling strength and r is the screening 
radius. The function $(^1,^2) vanishes for |£i,2| > 1, and 
is normalized by condition $(0,0) = 1. For simplicity it 
is taken to be an even function of its two variables. In 
terms of boson field operators 9i, the interaction acquires 
the following form, 

L/2 

TT Tr ^ f , , (x\-ri\a n 2 a~x 2 \ 
H int = V I dx 1 dx 2 $ [ , I > 



ro 



ro 



"-L/2 



d Xl 0i{x u n 2 a)d X2 9 2 (n 2 a,x 2 ) 



(3) 



Anticipating a Fourier expansion we introduce a a 2D 
basis of periodic Bloch functions which is constructed in 
terms of ID Bloch functions for non-interacting QB, 



where 



ip p ,q(x) = -j=e lqx f p (q,x), 



(4) 



and 



L/2 

H 2 = yY1 / dx 2 {g^l (n iai x 2 ) 

Ml -L/2 

- {d X2 2 (nia,x 2 )) 2 
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Here m enumerates the wires in the i-th array, Xi is a 
ID contiuous coordinate along e*, while (9i,iTi) are the 
conventional canonically conjugate boson fields (see, e.g., 
Rcf. [|l|]). The Fermi velocities vf 12 = v and the LL pa- 
rameters g\ 2 — 9 are taken to be the same for both 
arrays. The modulus of each component of a quasi- 
momentum in the first BZ does not exceed Q/2 where 
Q = 2n/a. Generalization to the case of different param- 
eters Vi,gi,di is straightforward. 

The interwire interaction results from a short-range con- 
tact capacitive coupling in the crosses of the bar, 

L/2 

Hint = ^2 J dx\dx 2 V{x\ — n\a, n 2 a — x 2 ) x 



Pi(xi,n 2 a)p 2 (nia,X2) 



(2) 



Here pi(r) are density operators, and V(ri— r 2 ) is a short- 
range interwire interaction. Physically, it represents a 
highly screened Coulomb electron-electron interaction 
between two charges located at the points ri = (x\, n 2 a) 
and r 2 = (n\a, x 2 ). In what follows we use an interaction 
of the form, 



f p (q,x) = exp ji sign (q) (-l) p+1 | Qxj . 

Here p, p' = 1,2,..., are the band numbers in a re- 
duced BZ, and q = {91,92} is the crystal quasimomen- 
tum, (\qi\ < Q/2). In momentum representation the full 
Hamiltonian (fy) then acquires the form, 
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with matrix elements for interwire coupling given by. 



(5) 
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(-i) p+1 h\) 



(6) 



are eigenfrequencies of the "unperturbed" ID mode per- 
taining to an array i, band p and quasimomentum q. The 
coefficients 



aipvp^ = sign(g 1(?2 )(-l)P+P'^o$ 



ipi 1 p' q? 



(7) 



are proportional to the dimensionless Fourier component 
of the interaction strengths 
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$ Wq = /da^2$(6,6)e-"-° (9l « 1+92 « 2) x 

x /p(gi,n)Ci)/p'(?2, 7-062 



< I ) 2p'lpq- (8) 



The Hamiltonian (||) describes a system of coupled har- 
monic oscillators, which can be exactly diagonalized with 
the help of a certain canonical linear transformation (note 
that due to spatial periodicity, it is already diagonal with 
respect to the quasimomentum q). The diagonalization 
procedure is, nevertheless, rather cumbersome due to the 
mixing of states belonging to different bands and arrays. 
However, it is seen from Eq.(0) that in the case ro <C a 
the dimensionless interaction a becomes effectively weak 
(numerical estimates are given below) and a perturbation 
approach is applicable. In this limit, the systematics of 
unperturbed levels and states is grossely conserved, at 
least in the low energy region corresponding to the first 
few bands. Indeed, as it follows from the unperturbed 
dispersion law (^J), the inter-band mixing is significant 
only along the high symmetry directions in the first BZ. 
The effect of interband interactions can be accounted for 
perturbatively in the rest of the BZ. Moreover, the mix- 
ing between modes within the same energy band is strong 
for waves with quasimomenta close to the diagonal of the 
first BZ. Away from the diagonal, this mixing effect can 
also be calculated perturbatively. 

In second order of perturbation theory the above men- 
tioned canonical transformation results in the following 
renormalizcd field operators and the corresponding renor- 
malized eigenfrequencies for the first array: 



Eaip2p'q^lpq^2p'q Q 
—2 ~72 y 2p'q, 

p, W lpq " W 2p'q 
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where 
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Corresponding formulas for the second array are obtained 
by replacing lp — ► 2p, and lp' — ► 2p' . 

For quasimomenta lying off the diagonal of the first BZ, 
all terms in these equations are non-singular. There- 
fore, the applicability of perturbation theory is related 
to the convergence of the series on the right hand side of 
Eq.©. Away from the first BZ boundary (|q| < Q/2) 
the following estimation is valid, 



*pq 



1 + I - 



V(b/2]Q)V 
(MO = CM), (2,1); p>l. 



Therefore, the correction can be approximately estimated 



as w^qS'q with 



S n = ^ a 2 



( Vogr 2 N " 
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For a short-range interaction, ro <C a, &ipi'p>q is a smooth 
function of p and p' . Therefore, the sum over p in Eq.(|ll|) 
can be replaced by an integral over the extended BZ with 
wave vector k whose components are 

fc l = (Zj +sign ( qi ){-l) p ' +1 [ Pl /2]Q. 

For |q| <§C Q one gets Sq = Sq(1 + o(l)) where, 



and 



-ir (feiCi+fc2?2) 



Finally, one arrives at the estimate, 



So 



Vbgjo 
fiv 



(12) 



where 



In the case of QB formed by nanotubcs, Vq is the 
Coulomb interaction screened at a distance of the or- 
der of the nanotube radiusE-tl i?o, so that Vq « e 2 /Ro and 
ro ~ Ro- Therefore, 



So 



hv 



(13) 



For carbon nanotubes one has0 v w 8 • 10 7 cm/sec and 
5 « 1/3. Then, $ can be calculated e.g. for a Gaussian 
model, 

*(£i,6) oc e ^ ( «i+«^ 2 « 1 « 2 cosv) . 

Taking for example </? = 7r/3, one finds that interaction is 
effectively weak if a ^> 4i?o- Therefore even slightly rar- 
efied QB already satisfies the desired condition. Thus, 
the dimensionless coupling constant a (0) is a small pa- 
rameter of the theory. Further calculations are carried 
out within an accuracy o(a 2 ). Yet, in the equations pre- 
sented below, only the first non-vanishing representative 
terms are retained (just for the sake of brevity). 
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Consider now modes with quasi-momenta near the diag- 
onal of the first Brillouin zone, but away from its bound- 
ary. In this case, the frequencies of the modes belonging 



to the same band coincide, wi pq 
the modes are strongly mixed: 



: ^2pq : 



J P<i 



Therefore 



?2pq 



'ipq 



(14) 



The corresponding eigenfrequencies are shifted from their 
bare values already in the first order in a 



~ 2 

W lpq 

~ 2 

W 2pq 
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(15) 
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FIG. 1. Solid and dashed lines of equal energy for interact- 
ing and noninteracting systems respectively. 

These results show that the quantum states of the 2D 
quantum bar conserve the quasi ID character of the 
Luttinger— like liquid in the major part of momentum 
space, and that the 2D effects can be calculated within 
the framework of perturbation theory. However, it fol- 
lows from Eqs.(Cl) that the bosons with quasimomenta 
close to the diagonal of the first BZ are the strongly mixed 
bare ID bosons. These excitations are essentially two- 
dimensional, and therefore the lines of equal energy in 
this part of BZ are modified by the 2D interaction (|l|). It 
is clearly seen that deviations from linearity occur only 
in a small part of the BZ. However, due to the absence of 
charge transfer, the Fermi surface of QB is still a "cross" 
obtained by superposition of two mutually perpendicular 
stripes of width \2kp\ around the lines k x = 0, k y = (cf. 
Rcf. HJ). The crossover from LL to FL behavior around 
isolated points of the BZ due to a single-particle hy- 
bridization (tunneling) for Fermi excitations was noticed 



in Refs. ||,[l(], where a mesh of horizontal and vertical 
stripes in superconducting cuprates was studied. 



III. CORRELATIONS AND OBSERVABLES 

The structure of the energy spectrum analyzed above 
predetermines optical and transport properties of the 
QB. Let us consider an optical conductivity aw (q, uj) 
whose spectral properties are given by a current-current 
correlator 




iiiq(t),j| lq (Q) 



(16) 



Here ji P q = \/2vgiri pci is a current operator. For sim- 
plicity we restrict ourselves to the first band. For non- 
interacting wires, the current-current correlator is re- 
duced to the conventional LL expressionlij, 

( iilq(*)iil'lq(°) ) = -Siws^iqSiri^iqO' 5 "' (17) 

with a metallic Drude peak 

au>(q,LJ > 0) = nvg6(u — u)iiq)6ii> . 

For interacting wires, where cti P i> p >q ^ 0, the correla- 
tors may be easily calculated after diagonalization of the 
Hamiltonian (|5|) by the transformations (Q) for q off the 
diagonal of the first BZ, or by the transformation (|l4| ) 
for q lying on the diagonal of the BZ. 

Consider first the optical conductivity for q far from the 
diagonal of the first BZ. In this case, the transformations 
for the field momenta can be obtained in a similar manner 
to the transformations @ for the field coordinates. As a 
result, one has: 



jllq(*),iil q (°) 

- 2ivg (1 - /3 q ) £ Uq sin(wn q i) - 
- 2ivgPqU> 2p < i sin(w 2 pqi), 

([jllq(*).illq(0)])« -2 i «5/5qX 

x (cjiiqsin(o)iiqt) - w 2 iqsin(a>2i q t)) 
where a q = an2i q and 

0!qWllq!jJ21q 



(18) 



/9q = 



ui 



llq 



UJ. 



21q 



Then one obtains 



(Tn (q, uj > 0) w Tivg (l-0%)6(u- cZ>n q ) + 

+ KVgfilfi (UJ - d> 2 pq) , 



(19) 
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<7i2 (q, oj > 0) w 7rv5/3 q x 

X (5 (CJ - 0>ll q ) - (5 (W - 0>21q) 



(20) 



The longitudinal optical conductivity ([ll]) (i.e. the con- 
ductivity within a given set of wires) has its main peak 
at the frequency wu q « v\qi\, corresponding to the first 
band of the pertinent array, and an additional weak peak 
at the frequency cj 2 i q ~ v\q2\i corresponding to the first 
band of a complementary array. It contains also a set of 
weak peaks at frequencies (I> 2 pq « [p/2]vQ (p = 2,3, . . .), 
omitted in Eq.([l9j) and corresponding to the contribu- 
tion from higher bands of the complementary array. At 
the same time, a second observable becomes relevant, 
namely, the transverse optical conductivity (j2(]). It is 
proportional to the interaction strength and has two 
peaks at frequencies <Dii q and ti>2iq in the first bands of 
both sets of wires. For |q| — * 0, Eq.(|l§|) reduces to that 
for an array of noninteracting wires, and the transverse 
optical conductivity vanishes. 

In case where the quasimomenta q belong to the diagonal 
of the first BZ, the transformations for the field momenta 
are similar in form to Eqs. (|l4|). The current-current 
correlation functions have the form 



WO,Jii q (o) 

- ivg pn q sin(wii q t) + <I> 2 i q sm(w 2 i q t)] , 

([ilq(f),i|q(0)]) « 

- ivg (<Dii q sin(cjn q <) - tl> 2 i q sin(ai 2 iqt)) . 
and the optical conductivity is estimated as, 



(21) 



CTn(q,w > 0) 



cri2 (q, w > 0) 



irvg 



6(U- Wnq) +5((JJ- UJ 2 lq) 



(22) 



5(w- wn q ) -5(u- w 2 i q ) . (23) 



Here, the index i of the array is omitted because the 
frequencies for both arrays coincide, cji pq = w 2pq - The 
longitudinal optical conductivity (^2|) has a split dou- 
ble peak at frequencies tDii q and Waiq, instead of a sin- 
gle peak. Again, a series of weak peaks (omitted in the 
r.h.s. of Eq.(|22])) occurs at frequencies w pq correspond- 
ing to contribution from higher bands p = 2, 3, 4, The 

transverse optical conductivity (p3|), similarly to the non- 
diagonal case ( po| ) , has a split double peak at frequencies 
a>ii q and w 2 iq- 

One of the main effects specific for a QB is the appear- 
ance of non-zero transverse momentum-momentum cor- 
relation function. In space-time coordinates (x, t) its rep- 
resentation reads, 



This function describes the momentum response at the 
point (0, X2) of the second array at the moment t caused 
by initial (t = 0) perturbation at the point (xiO) of the 
first array. Standard calculations similar to those de- 
scribed above, lead to the following expression, 



Gi 2 (x;i) 



Vor 2 
2irhva 



dkidki'&kxkzkikz x 



x sin(fcixi) sin(fc2£2 



&2 sin(fc2i>i) — fci sin(fciut) 
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FIG. 2. The transverse correlation function Gi2{xi, X2;t) 
for ro = 1 and vt — 10. 

This correlator is shown in Fig [2] for xi$ > 0. Here the 
non-zero response corresponds to the line determined by 
the obvious kinematic condition x\ + x 2 = vt. The finitc- 
ness of the interaction radius slightly spreads this peak 
and changes its profile. 

Further manifestation of the 2D character of QB system 
is related to a possible periodic energy transfer between 
the two arrays of wires. Consider an initial perturba- 
tion which, in the system of non-interacting arrays, cor- 
responds to a plane wave propagating within the first 
array along the ei direction, 

vpo 



(6i(xi,n 2 a;t)) - _ 

V 2wii q 

(9 2 (n 1 a 1 x 2 ;t)) = 0, 



sin(gia;i + q 2 n 2 a - wn q t), 

(25) 



G 12 {x,t) = ([^1 (xi, 0; 0,^2(0, £25 0)]). 



(24) 



(po is the charge density amplitude). If the wave vector q, 
satisfying the condition |q| << Q/2, is not close to the 
first BZ diagonal, weak interwire interaction a slightly 
changes the (9±) component and leads to the appearance 
of a small (9 2 ) ~ a component. But for q lying on the di- 
agonal, both components within the main approximation 
have the same order of magnitude 
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h {xi,n 2 a;t) 



vpo (I , 

COS -ftqWlqi 



V2lu 



lq 



x sin(gia;i + (72^20 — wi q i), 



(26) 



vpo 



V2lu 



Sill I -OlqUJlqt 



lq 



x cos(ginia + 172X2 — ^iqt). 

This corresponds to a 2D propagation of a plane wave 
with wave vector q, modulated by a "slow" frequency 
~ auj. As a result, an energy is periodically trans- 
fered from one array to another during a long period 
T ~ (aui)~ 1 . These peculiar "Rabi oscillations" may be 
considered as one of the fingerprints of the physics ex- 
posed in QB systems. 
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IV. CONCLUSION 



In conclusion, we have demonstrated that the energy 
spectrum of QB shows the characteristic properties of 
LL at \q\,u — > 0, but at finite q the density and momen- 
tum waves may have either ID or 2D character depend- 
ing on the direction of the wave vector. Due to inter- 
wire interaction, unperturbed states, propagating along 
the two arrays are always mixed, and transverse compo- 
nents of correlation functions do not vanish. For quasi- 
momentum lying on the diagonal of the Brillouin zone, 
such a mixing is strong and transverse correlators have 
the same order of magnitude as the longitudinal ones. 
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